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ABSTRACT 

We propose a method of distance determination based on the internal structure and 
dynamics of disk galaxies. The method relies on the universal luminosity profile of a stellar 
disk represented by an exponential law. Calibrating nearby galaxies with known distances, 
it is found that the scale length of the disk is tightly correlated with the specific combination 
of central surface brightness and rotational velocity at a characteristic radius of 2.15 scale 
lengths from the center. This suggests that the scale length of the disk may be used as an 
indicator for extragalactic distance scale. The application of this relation to M51 and MlOO 
allows us to arrive at the distances of about 6 Mpc and 14 Mpc, respectively, implying a 
Hubble constant of i/o = 92 ~ 94 km s"^ Mpc"^ 

Subject headings: distance scale — galaxies: kinematics and dynamics — photometry 
— spiral — structure 

1. INTRODUCTION 

Accurate measurement of the distances to galaxies is fundamental for understanding 
how fast the universe expands as well as how much the universe is aged. For distant galaxies 
beyond a few megaparsecs, the reliable distance estimates from Cepheids are unavailable, 
and it is customary to construct a cosmic distance ladder using distance-free empirical 
correlations among the structural and photometric properties of galaxies (e.g. Tully & 
Fisher 1977; Faber & Jackson 1976; Tonry & Schneider 1988; Pierce & Tully 1988; van den 
Bergh 1992). 

Most of empirical correlations for galaxies originate from the virial theorem which 
yields a three-parameter relation of (luminosity) oc (velocity)^ x (length) if the mass-to- 
light ratio is constant. For some particular systems like spiral galaxies, the virial theo- 
rem possibly degenerates into two-parameter relations, because the central surface bright- 
ness of stellar disk is known to be in a narrow range, yielding an additional constraint 
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{luminosity) I {I engthY ^ constant (Freeman 1970). The above dimensional argument ex- 
plains the background for any souglit-after relation for galaxies. It is the tightness of the 
relation that provides a reliable distance indicator, and the tightness itself depends upon 
how well the correlated quantities can be defined among many different forms which have 
the same dimensionality. 

A widely used method of measuring the distances to spiral galaxies by TuUy & Fisher 
(1977) is based on either the luminosity - velocity relation or the velocity - diameter relation 
(the former is more popular and is called the TuUy-Fisher relation). More specifically, 
TuUy & Fisher used the total luminosities, HI velocity widths, and apparent diameters. 
We note that the total luminosity represents the luminous matter in both the bulge and 
disk. On the other hand, the velocity width, which corresponds to the maximum rotational 
velocity at larger radii, is more closely related to a dark matter in halo and can be affected 
by distortion due to a companion galaxy. Thus, the correlation was made between the 
quantities which represent the different components. When the diameter of a galaxy at a 
fixed isophote is used as a dimension of length, brighter galaxies are biased in favor of larger 
intrinsic diameters. All these factors increase the dispersion of the TuUy- Fisher relation, 
which enhances errors in the distances to be estimated. 

One of the conditions for establishing an ideal distance indicator is that the correlated 
quantities are chosen on a well-defined physical basis, so that some uncertainties associated 
with the correlation can be controlled (Aaronson & Mould 1986). We aim to find a more 
precise distance indicator by using a sensible combination of observable quantities which 
represent only the disk with little contamination from the bulge and halo. 

In §2 we describe the basic scheme of our method. In §3, using the observed quantities 
for local calibrator galaxies, we derive the fundamental relation which can be used as 
a distance indicator. The advantages of our method are discussed in §4. Applying our 
method to the galaxies of M51 and MlOO, we estimate their distances and a corresponding 
Hubble constant in §5. 

2. METHOD 

The radial luminosity profile of a galactic disk is well fitted by an exponential function, 
/(r) = Iq exp( — r//i), where Iq is the extrapolated central surface brightness and h the scale 
length of the disk (Freeman 1970). Under the assumptions that the disk is in centrifugal 
equilibrium and the mass-to-light ratio {M j L)d of the disk is constant, such an exponential 
profile gives a rotation curve that attains a maximum Vd/ j2'KGhIo{M/ L)d — 0.62(= /) 
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at r jh — 2.15. In an actual galaxy, the overall shape of its rotation curve would differ 
from what is expected solely from the exponential disk unless a bulge and/or dark halo 
constitutes a minor component. However, such deviation would be rather small aX r/h — 
2.15 where / is close to 0.62, because in most cases this radius is well outside the bulge and 
well inside the dark halo. Thus, the rotational velocity observed at r/h — 2.15 manifests 
the dynamical state of a luminous matter of the disk. 

Given the distance-independent observables /q and Vj. for a galaxy, we can derive the 
intrinsic scale length h using the formula which should hold for the disk: 



Then, comparing this intrinsic length h to the observed length h measured in angular units, 
it is straightforward to determine the distance to the galaxy under consideration. 

In the above formula, the quantity s(M, L) contains the mass-to-light ratio Mj L oc L'^ 
which would vary from galaxy to galaxy with some power of its luminosity (e.g., Salucci, 
Ashman, & Persic 1991). It is difficult to predict the explicit form of s(M, L) from purely 
theoretical grounds: we shall assume that s is represented as a functional of the distance- 
independent variable VJ/Iq. In the next section we constrain a possible form of s using 
the local calibrator galaxies whose distances have been estimated from Cepheids and other 
methods. 



Table 1 tabulates the adopted data for the local calibrator galaxies. The distance 
moduli (DM) (Col.(3)) to the galaxies, which have been estimated directly from Cepheids, 
are taken from de Vaucouleurs (1993). For other galaxies the DMs are the average of a 
number of independent DM determinations. DMs based on the brightest stars are taken 
from Tammann (1987). The DM of NGC55 is from the indirect distance estimate using 
carbon stars (van den Bergh 1992), and the DM of NGC5585 is assumed to be equal to 
that of MlOl because they belong to the same group. For M81, we adopt the long-distance 
scale (Tammann 1987), but this does not change the following result significantly (see 
below). Columns 4 and 5 contain the disk scale length h (arcmin) and the central surface 
brightness /ig (mag arcsec"'^) in the B band, together with their la errors. Since the 
decomposition of luminosity profile into the exponential disk includes some uncertainty, 
we have re-determined h and /io as well as their errors. The references given in Col. (6) 
originally published the luminosity profile. Presented in Cols.(4)-(5) are the average values 
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of our own estimates with those from the references. Column 7 presents tlie observed 
line-of-siglit velocity at r — 2.15/i, while Column 8 for the inclination i which has been 
evaluated by analyzing the velocity field as given in the references (Col.(9)). 

The face-on central brightness /y {Lbq pc^^) is derived from /io, after correction for 
the absorption and inclination. The absolute i?-magnitude of the Sun is adopted as 
Mbq = 5.48. The corrections for the Galactic and internal absorption are taken from Third 
Reference Catalogue of Bright Galaxies (RC3, de Vaucouleurs et al. 1991). We have also de- 
rived the inclination of the optical disk from the standard formula, cos'^ i — (q^ — q^) l{l—q^), 
where q is an axial ratio at an isophote of 25 mag arcsec"^ (RC3) and q^ is fixed to be 0.2. 
This photometry-based inclination is not necessarily equal to the velocity-based inclination 
given in Table 1. To keep consistency in the analysis below, we use the photometry-based 
inclination in correcting both Vd and /q. 

1 /2 

Figure 1 shows the data in the (log /i, log V^^/Zq ) diagram. Filled circles represent the 
galaxies for which the distances are estimated directly by using Cepheids, while empty cir- 
cles for those estimated by other means. The tight correlation between these two quantities 
is apparent, and this correlation is fitted by a linear relation: 

logh^alogVdlll'' + h , (2) 

where a and h are constants to be determined by the least-squares fit to the data points. 
We use only the data with direct distance estimates (filled circles) and weight them by an 

1 /2 

inverse of the sum of errors in log h and log Vdj Iq ■ The best fit line is shown by solid line in 
Figure 2. The slope of this line is a = 1.045 ±0.072 and its intercept is 5 = -0.658 ±0.015. 
We call this a standard Case A in Table 2. 

There is a freedom for choice of some controversial distance estimates (short or long 
scale) as well as a freedom for use of an alternative scheme of absorption correction based on 
A Revised Shapley-Ames Catalog of Bright Calaxies (RSA, Sandage & Tammann 1981). In 
order to examine the sensitivity of the least-squares fit to these modifications, we consider 
five different Cases B to F and their results are tabulated in Table 2. The label "T" 
in the column DM of Table 2 stands for exclusive use of Tammann's (1987) distance 
estimates except for NGC3109 (Sandage & Carlson 1988; Capaccioli et al. 1992) and 
NGC55 (van den Bergh 1992). The label "vdB" is for van den Bergh's (1992) distance 
estimates except for NGC247, NGC253 and NGC7793 (Tammann 1987) and NGC3109 
(Sandage & Carlson 1988; Capaccioli et al. 1992). When compared to Case A, the most 
noticeable change appears in the scale lengths of M81 and MlOl, as shown by triangles in 
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the inset of Figure 1. 

The values of a and b in Cases B to F are well within la error of the corresponding 
values in Case A. The reason is that a less weight is assigned to these big galaxies because 
visible fine structures such as spiral arms affect reliable determination of the background 
exponential disk. In other words, the galaxies with regular luminosity profiles are obviously 
given a high weight in the least-squares analysis. We further note that the fit is also 
insensitive to the absorption-correction scheme, either RC3 or RSA. 

There should be another correlation between h and Vd if /o resides in a narrow range 
and is approximated as a constant. We have confirmed this correlation but with a much 
larger dispersion. As shown in Figure 1, the fact that inclusion of /q as a third quantity 
greatly improves the correlation justifies the specific way of combining Iq with Vd in Eq.(l). 

4. DISCUSSION 

4.1. IncUnation Correction 

One of the remarkable points is that the nearly face-on galaxy MlOl with i — 27° is 
close to the line of the best fit in Figure 1, which suggests that our method is also applicable 
to galaxies with such small inclinations. This is contrasted with the TuUy-Fisher relation 
for which the sample is restricted to highly-inclined galaxies with i > 45° to minimize 
the line-of-sight velocity correction, but with i < 85° to avoid large correction for internal 
absorption. 

Our method uses the relation of the inclination-independent h against the variable 

1 /2 

Vdj Iq which depends on i as follows: 

log % - log (t^)" - \ log{(sinzr(coszr--(T)| ^ (3) 

where the superscript uncorr denotes the quantity without inclination correction. The 
factor a(T) stems from the internal absorption, and it depends on galaxy type index T. 
The explicit form of a(T) is given in RC3. Equation (3) is used to figure out how the error 

1 /2 

in i propagates to logVd/lQ , Given Ai — 5°, for example, and oc{T) — 1.4 for an average 

1 /2 

type index T = 6.5 of local calibrators, we calculate the propagated error AlogVd/ Iq 
and the result is shown by the solid line in Figure 2. 

1 /2 

The above exercise shows that AlogVd/ 1"/ for MlOl {i ~ 27°) is as small as ±0.06 
and it is equal to that for highly- inclined galaxies with i ~ 85°. Remembering that the 
slope in Eq.(2) is a ~ 1 as shown in Table 2, we can evaluate the error in the resulting 

1/2 

determination of a Hubble constant Hq using AlogVd/ Iq' ~ Alog/i ~ Alogifo. When 
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the sample is restricted to 45° < i < 80° as in the Tully-Fisher method, it turns out 
that the uncertainty of log Hq due to the inclination correction is well within ±0.02 and is 
practically negligible. 

4.2. The Mass-to-Light Ratio of the Disk 

The tight relation (Eq.(2)) which we found, guided by a theoretical argument, relies 
heavily on the local calibrator galaxies. Therefore, the derived relation still needs a physical 
explanation. The dynamical equation in Eq.(l) is exact for a disk. The scaling relation 
h oc {Vd/Io^'^Y (Eq.(2)) should be equivalent to f{M/L)D oc {Vd/ lo^^'^-^'^ ■ It is this 
regularity of the internal structure and dynamics of the disk that exists behind the tightness 
of the derived relation. Our result of a ~ 1 as well as / ~ 0.6 indicates that the mass-to- 
light ratio of the disk varies from galaxy to galaxy as (M/L)^) oc Vd/Io^"^ oc {Ljj/ IqY^'^. It 
is interesting that the value of a is almost exactly equal to an integer of unity (see Table 
2). If confirmed from more extensive data, it is tempting to explain its origin within the 
context of galaxy formation and evolution. 

4.3. Summary of the Proposed Method 

Our method uses the three quantities of the rotation velocity, the scale length of the 
disk, and the central surface brightness of the disk which are designed to represent the disk 
with little contamination of the bulge and halo. 

The rotational velocity Vd at a relatively inner radius of r/h — 2.15 is contributed 
maximally by the exponential mass distribution of the disk. The velocity measured at this 
radius is almost irrelevant to how much the outer velocity field is distorted due to tidal 
effects. The scale length h and the central surface brightness /io are determined through 
the decomposition of the luminosity profile, in a coupled way that the larger h results in 
fainter /Iq. These quantities are basically independent of the outer cut-off of the disk which 
results from the isophotal threshold above which the flux is measured. 

Local illumination of young stars in spiral arms affects the proflle decomposition, but 
the resulting h and /io change their location, more or less, along the correlation in Fig.l. 
We see this from the inset where the results of two different decompositions for M31 are 
shown. It is evident from this flgure that the uncertainty of the proflle decomposition 
hardly violates the correlation derived in §3. 

On grounds of dimensional argument, our method is transformed to the virial theorem. 
The application to distance estimates is traced back to Oepik (1922). The present paper 
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basically demonstrates the virial method, but from a practical viewpoint of utilizing modern 
data to determine distances with precision. The merits of our method are: (1) the result is 
very insensitive to the inclination correction, so the restriction on i is greatly relaxed, and 
(2) the correlated quantities have a physical basis and are well defined and unambiguous. 
The apparent drawback is that detailed observations of photometry and kinematics are 
required. However, observation with modern equipment can provide the required data for 
the precise distance measurements advocated here. 

5. APPLICATION OF THE METHOD 
The relation we propose to determine the distance is (Case A): 

log h = 1.045(±0.072) log Vd/Io^^ - 0.658(±0.015) (x^ = 0.0062) , (4) 

where the units are /i(kpc), V^(km/s), and Iq{Lbq pc^^). We apply this relation to the 
Whirlpool galaxy, M51, and a member of the Virgo Cluster, MIOO. The basic quantities 
of M51 are h — 111.60 arcsec, /io = 21.565 mag arcsec"^, Vd — 164.23 km and 
i = 20° (Schweizer 1976; Burkhead 1978; Elmegreen & Elmegreen 1984; TuUy 1974), and 
those for MIOO are h — 56.60 arcsec, /io = 21.335 mag arcsec^^, Vd — 213.80 km s^-^ 
and i — 27° (Schweizer 1976; Benedict 1976; Elmegreen & Elmegreen 1984; Kodaira et 
al. 1986; Distefano et al. 1990). The RC3 scheme has been used for the correction of 
Galactic and internal absorption, and the inclination has been estimated from the analysis 
of velocity field of the galaxy. Note that the values of /io and Vd are those corrected for the 
absorption and inclination. 

For M51, the velocity-based inclination i — 20° differs from the photometry-based incli- 
nation taken from RC3, most probably due to the strong interaction with the companion 
galaxy. Since the absorption correction is negligible in a nearly face-on galaxy, we have 
adopted i — 20°, arriving at the distance of 5.97i;'j;g8 Mpc. The recession velocity of M51 
is adopted as 551 km s~-^ from RC3, then we obtain the Hubble constant Hq — 92.31^15^1 
km s^-^ Mpc^-^. 

For MIOO, the velocity-based inclination i — 27° results in the distance of 14.25^i;g7 
Mpc, while 13.88ii;g3 Mpc for the photometry-based inclination, thereby arriving at the 
mean distance of lAMt\f^ Mpc. The systematic recession velocity of the Virgo Cluster is 
estimated as 1316±120 km s^-^, after corrected for the Virgocentric infall taken from Pierce 
& TuUy (1988). With this recession velocity we obtain Hq = 'dZ.Qt^^^ km s"^ Mpc'^ 
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We are grateful to the referee for his valuable comments on our method, and to B.A. 
Peterson and N. Visvanathan for comments on the manuscript. 
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TABLE 1 
Local Calibration Data 
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"Cepheid distance estimates are available. 

Notes.— Col.(3):Distance moduli from de Vaucouleurs 1993, except for NGC247, NGC253, NGC7793, 

M81 (Tammann 1987) and NGC55 (van den Bergh 1992). 
Col.(4):Disk scale lengths. 

Col.(5):Extrapolated central surface brightnesses of the disks in the B band (uncorrected for inclination 
and absorption). 

Col.(6):References for surface photometry parameters: (l)Elson et al. 1987; (2)Kent 1987; 

(3) de Vaucouleurs 1958; (4)de Vaucouleurs 1959; (5)Carignan 1985a; (6)Carignan 1985b; (7)Puche et 
al. 1991a; (8)Puche et al. 1991b; (9)Pence 1980; (lO)Okamura et al. 1977; (ll)Schweizer 1976; 
(12)Yoshizawa & Wakamatsu 1975; (13)Cote et al. 1991; (14)Okamura et al. 1976. 

Col.(7):Line-of-sight velocities at a radius oi r/h — 2.15 (uncorrected for inclination). 
Col. (8) inclinations estimated from analysis of velocity fields. 

Col.(9):References for rotation velocities:(l)Elson et al. 1987; (2)Emerson 1976; (3)Rogstad et al. 1976; 

(4) Carignan 1985a; (5)Kent 1987; (6)Puche & Carignan 1991; (7)Rogstad et al. 1979; (8)Carignan & 
Puche 1990b; (9)Carignan & Puche 1990a; (lO)Puche et al. 1991a; (ll)Puche et al. 1991b; 
(12)Shostak 1973; (13)Rots 1975; (14)Visser 1980; (15)Cote et al. 1991; (16)Comte et al. 1979. 
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TABLE 2 

Least-Square Fitting Parameters 



Case 


DM" 


Absorption'' 


slope 


intercept 
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A 
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1.045±0.072 


-0.658±0.015 


0.0062 
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-0.668±0.028 


0.0173 
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T 
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1.031±0.077 
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0.005C 
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T 
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1.041±0.136 


-0.641±0.028 


0.0148 


E 


vdB 


RC3 


1.011±0.118 


-0.624±0.026 


0.0089 


F 


vdB 


RSA 


1.022±0.151 


-0.630±0.031 


0.0018 



"Distance moduli mainly from de Vaucouleurs 1993 (deV), Tammann 1987 (T), 
and van den Bergli 1992 (vdB). For more explanation see text. 
''Absorption corrections from RC3 or RSA catalogue. 
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Figure Captions 

FIG. 1 Log-Log plot of the scale length of a disk h (kpc) against tlie combination of 
the rotational velocity (km s^-^) at r/h — 2.15 and the extrapolated central surface 
brightness /y {Lbq pc^^) for Case A. Filled circles denote the local calibrator galaxies 
for which the direct distance measurements are available, while empty circles denote 
the calibrators with only indirect distance measurements. The solid straight line is the 
result of weighted least-squares fit to the filled circles. The tnset: Triangles show how 
the locations of M81 and MlOl in the diagram change for Case E with different distance 
moduli, compared to those in Case A (filled circles). Empty circle for M31 corresponds 
to one of the different decompositions for its photometric data from which the mean and 
error are estimated (filled circle). Each line is the result of weighted least-squares fit for 
Case A (solid), B (dotted), C(short-dashed) and E(long-dashed). Note the very small 
difference among the fitted lines. 

1 /2 

FIG. 2 The error in logy^j/Jg arising from the uncertainty of inclination Ai — 5°. The 
dependence of internal absorption correction on axial ratio is written as alog(major 
axis/minor axis) with a — lA for the average type index T = 6.5 for the local calibrator 
galaxies in Table 1. 
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